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We put forward an interpretation of scalar quantum field theory as relativistic quantum mechanics 
by curing well known problems related to locality. A probabilistic interpretation of quantum field 
theory similar to quantum mechanics is difficult if particle localization is defined using the Newton- 
Wigner position operator as it is non-local and non-covariant. An alternative bilinear covariant 
position operator is discussed which incorporates a time operator that can be exponentiated to 
a unitary operator. Moreover, it satisfies an algebra that unifies special relativity and quantum 
mechanics and has the same form for particles with spin. Higher power position operators are 
derived which yield Heisenberg's uncertainty relations. Our ideas are illustrated with a relativistic 
wave function whose probability density can be perfectly localized. 



I. OUTLINE 

In this paper we present a consistent first quantization of scalar quantum field theory. The term "first quantization" 
refers to the formulation of relativistic quantum field theory as relativistic quantum mechanics with a Fock space of 
single particle wave functions that have a Copenhagen-like probabilistic interpretation. In appendix \^ we briefly 
address the standard relativistic quantum mechanics interpretation of quantum field theory and its problems. We 
recall that in standard textbook treatments on scattering theory single particle wave functions are identified with 
positive frequency solutions of the Klein Gordon equation. The restriction to positive frequency solutions implies that 
the time evolution is governed by the Feynman propagator. Even though the wave function interpretation of positive 
frequency solutions is used to compute scattering probabilities and cross sections this interpretation is generally not 
deemed fundamental as it suffers from two well known problems related to non-locality, 

1. Time evolution as described by the Feynman propagator violates Einstein causality. 
In other words, it causes wave functions to spread outside of the lightcone. 

2. The conventional Newton- Wigner position operator [Tj is inherently non-local and non-covariant. 

In appendix |X] we review the first problem and in appendix [B] we summarize recent work Q which resolves it. We 
show that adding negative frequency modes to the wave function restores Einstein causality, precisely as in classical 
field theory. As in [2| we emphasize that negative frequency modes actually have positive energy as dictated by the 
energy momentum tensor. The time ordered time evolution of wave functions which include negative frequency modes 
is governed by the half-advanced half-retarded Wheeler propagator and therefore replaces the Feynman propagator^. 
The inclusion of negative energy modes moreover implies that the Klein Gordon field itself is interpreted as a wave 
function. We do so by recalling a form of the probability density from \^ that was inspired by @ and 8^ which is 
constructed from the field and the Hadamard two point function^. 

Our emphasis on causal methods for quantum field theory is partially motivated by the pivotal role causality plays in 
a recent nonperturbative approach to quantum gravity called Causal Dynamical Triangulations (CDT). This method 
was introduced in 9] where two dimensional world sheets were quantized while preserving a causal foliation. This 
work has been generalized to 3 -I- 1 dimensions using computer simulations, for recent work see for instance [Tol. [Tl|. 
Additionally, the causal world sheet perspective was further developed to include interactions by topology changes 

MM- 

In section we introduce an, as far as we are know, new and very simple form of the probability density which is 
constructed from the field and its Hilbert transform, a transform that is widely used in signal analysis. We interpret 



^ Within our conventions the Wheeler propagator is the real part of the Feynman propagator and it is also known as the "principal part 

propagator" , see for instance 0| , due the appearance of the principal part in its Fourier representation . The Wheeler propagator 

can also be thought of as being half the time ordered commutator 2] . 
^ In the cosmological literature the Hadamard function is referred to as the statistical propagator, and in second quantization it is given 

by the anti-commutator. As observed in we emphasize however that the Hadamard function is not a propagator in the sense that 

it does not govern the dynamical propagation of a field in time, see sections IB 31 and HD of this paper. 
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this new incarnation of the probabiUty density as a representation in a different basis of the positive/negative frequency 
form imphed by the innerproduct presented by Woodard ^] and Halhwell- Ortiz Q. 

In section IIIII we address the second source of non-locaUty in standard quantum field theory which stems from 
the localization of relativistic particles through the Newton- Wigner position operator [H. We discuss an alternative 
bilinear position operator which shares many of the desirable properties of the Newton- Wigner operator but does not 
possess any non-local terms, both in momentum and position space. Furthermore, the bilinear position operator is part 
of a covariant space-time position operator, unlike the Newton- Wigner operator. Additionally the bilinear operator 
is not by itself related to dynamics and is therefore a generalization of the Schrodinger picture position operator 
in non-relativistic quantum mechanics. The Newton- Wigner operator on the other hand satisfies a Heisenberg like 
equation. 

A particularly distinguishing feature of the bilinear operator is that it generalizes to higher powers of the coordinates. 
Such higher power position operators which are Hermitian with respect to relativistically normalized wave functions 
have to the best of our knowledge not been derived before. They allow us to re-derive the Heisenberg uncertainty 
relations from quantum field theory viewed as relativistic quantum mechanics. It is quite remarkable that this corner 
stone of quantum mechanics is no longer present in quantum field theory if one insists on using the Newton- Wigner 
operator. Consequently, maiiy authors have abandoned a quantum mechanics interpretation of quantum field theory 
entirely, see for instance [1Jt13- We show that using the bilinear position operators on wave functions that include 
negative frequencies restores the quantum mechanical interpretation of free quantum field theory. It is expected that 
a similar picture to holds in interacting field theory provided the interaction can be treated perturbatively. 

To conclude our discussion on the position operator, the observation is made that unlike the Newton- Wigner 
operator the bilinear operator is purely a position operator in the sense that it does not contain any other generators 
of the Lorentz group. This has the important consequence that the form of the operator does not depend on the spin 
or mass of the fields under consideration. 

In section IIVI we discuss a real solution to the Klein Gordon equation that is normalized with respect to our 
probability density using the Hilbert transform. It behaves qualitatively similar to a gaussian wave function and can 
therefore be perfectly localized if the width approaches zero, where the width is defined using our bilinear position 
and position squared operators. Infinitely precise localization requires infinite energy however, which is consistent 
with physical expectations related to Heisenberg's uncertainty relations. We contrast our states with the "localized 
states" of Newton and Wigner [T| which are neither normalizable nor perfectly localizable. 



II. THE HILBERT TRANSFORM AND THE PROBABILITY DENSITY 



The position space form of the probability density that is used in textbook treatments of scattering probabilities 
can be written as follows'^, 

P{x) ^ cl)l{x)idoM^)- (1) 

We refer to this as the Klein Gordon probability density. If one adds negative frequency modes^ following Q and @ 
one can find the equivalent form, 

Pi^) = I (^c^*+{x)idoMx) - <P-{x)ido<j>-ix)^ . (2) 

As discussed in [3], the Hadamard function can be used rewrite this in a manifestly real probability density on the 
space of real solutions of the Klein Gordon equation in position space. 

P{x) = l- f d^x'cl>{x)dtAH{x,x')dt,cb{x') , (3) 
2 J t=t' 

The Hadamard form of the probability density ([3]) contains an integral however which makes it slightly unnatural as 
a Born-like probability. Therefore we present an, as far as we know, new form of the probability density. Similar to 



^ For a discussion of the probabilistic interpretation used in conventional scattering theory see appendix 1X1 
* For details on the addition of negative frequency modes see appendix IB] 
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the Hadamard form, it is defined in terms of real fields but does not contain an integral in its definition^, 

P{x)^^cP{x)do<PH{x), (4) 

where 4'h{x) denotes the Hilbert transform of in the time variable. The Hilbert transform can be written as a 
principal part integral 

oo 
— oo 

or more explicitly as 

oo 

0.^1imlL-'^(^ + ^'--^-^(^-^'-). (6) 

€ 

Note that with these definitions the Hilbert transform is performed in the time direction. One can nonetheless extend 
the Hilbert transform in a Lorentz covariant manner to any time like direction, see [3] . Lorentz invariance of the 
Hilbert transform can furthermore be inferred from the fact that it only affects the sign of the energy modes of the 
field, which inherently is a Lorentz invariant operation. 

The position space representation of the probability density (j4]) for real fields has a simple, manifestly real form. 
At first sight it seems surprising that we can define a probability density entirely with real variables. A closer look 
however, reveals that although the definition of the Hilbert transform ^ does not involve complex numbers it defines 
a complex structure on the space of solutions in the sense that if it is applied twice it gives back the negative of the 
original function, 

Jh'P = -0, Jh4> = 4>H- (7) 

The Hilbert transform is well known in signal analysis and one of it's essential qualities is that it fiips the sign of the 
negative frequency modes of a signal. 



, = 0+ + <?!)_ =5R(A+. (8) 



From this property one can easily see that our definition of the innerproduct is equivalent to the form ([2]). The Fourier 
decomposition of the field and its Hilbert transformed field can thus be written as, 

= / (03'^(^' + m')f{k)e^''^ <Ph{x) = I J + m2)e(fc0)/(fc)e^'=% (9) 

where /(fc) is the wave packet in momentum space and if we define the following 



M^) = J J^Sie + m^)/^(fc)e^^-- = J j^S{k' + m^)/^(fc)e-*^^ (10) 
we see that 

fnik) = *e(fc°)/(fc), /^(fc) = -ze(fcO)/*(fc). (11) 

Our conventions for the Hilbert transform are such that Jh cos{ijjt) — sin(a;i). We stress that the Hilbert transform 
of a real signal is again a real signal, which means that c/jh is real. This should be clear from the fact that it can be 
written as minus the imaginary part of the positive frequency field. The rationale is often actually reversed and the 



^ A similar form of the innerproduct is discussed in the context of PT-symmetric quantum mechanics by A. Mostafazadeh Instead of 
Hilbert's integral transform the author utilizes a non-local differential operator C which is only defined for free theories. The operator 
C can be related [6^ to the charge grading operator that appears in the inner product of PT-symmetric quantum mechanics introduced 
by Bender, Brody and Jones [gtII . Because of the similarities between the Hilbert transform and C one might therefore deduce that also 
the Hilbert transform can be given such an interpretation. 
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Hilbert transform is used to, uniquely, define the complex valued positive frequency signal from the real valued input. 
From this point of view the positive frequency component (x) can be considered as an analytic signal representation 
of the field 0. 

From the Hadamard form of the probability density ([3]) it is clear that the Hilbert transform of the field can be 
obtained by acting on it with the Hadamard propagator, 

<j)H{x)= f d^x'AH{x,x')dt'cf>{x'), (12) 



We gave here a presentation of our probability density (jl]) as a manifestly real form of ([2]) which is the negative 
frequency generalization of the Klein Gordon density One does not need to make such an interpretation related 
to Fourier modes however, one can just take our probability formula ^ at face value. It is possible to compute the 
time evolution of the field through the Wheeler propagator (jBl[) , compute the Hilbert transform of the time evolved 
field and insert both in our probability density In this way one can formulate quantum mechanics entirely without 
complex valued wave functions^ . We stress furthermore that from the same point of view the Klein Gordon equation 
for real scalar fields is the relativistic generalization of the Schrodinger equation. Moreover, since the time evolution 
is governed by the Wheeler propagator instead of the Feynman propagator the evolution is manifestly causal, see Q 
and appendix [Bl 

A. Particles and anti particles: complex scalar fields 

In the discussion above we focused on the real scalar field and how it can be interpreted as a real valued quantum 
mechanical wave function. To properly discuss particles and antiparticles however it is instructive to examine a 
complex scalar field. We do so by writing the complex field as a combination of two real fields, 

ip = (pi + i(f>2- (13) 

Within our formalism ip and ip* are the wave functions for a charged particle and an anti-particle respectively. The 
main physical difference with the real scalar field is that while a real field is neutral, a complex field carries electric 
charge and the charge density is given by 

p = <p*{x)idtp{x) = 029*01 - 0i9f02 = 2(/)29f(/>i, (14) 

where the diagonal terms have disappeared because the 0i and 02 fields are electrically neutral separately. The naive 
generalization of our probability density of real fields Q is not real, 

p*{x)dtp^H{x) = i (^0iat0f + 02 9*0^ + ^01 9t0f-z02 a t0f^ . (15) 

The proper probability density for complex scalar fields is just the real part of the above, 

p{x) = np*{x)dtPH{x) = i (^0iat0f + 02^*0^^ , (16) 

in this form the probability density is manifestly real as required for a probability density. If we alternatively write 
it terms of positive and negative frequency wave functions we see that it has the exact same form as for real fields, 

P{x) = \ (^^*+{x)idt^+{x) - p*_{x)idtp-{x)^ , (17) 
and is the generalization of the Klein Gordon probability density for complex fields, 

PKcix) =P'*+ix)idtp>+ix). (18) 

A consistency check that our probability density formula (jl6l) is indeed a probability density for a charged particle 
is that when replacing p with its charge conjugate p* the density is left invariant. In other words, the probability 



The wave functions are real for neutral scalar fields, charged scalar fields are modeled as complex fields. 
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is the same for positively and negatively charged particles as it should be. Similarly, acting with the same charge 
conjugation operation on our formula for the charge density for wavefunctions that include negative energy modes 
(fT4| has the physically required effect of changing its sign. 

In this subsection we have constructed a manifestly real probability density on the space of complex solutions to 
the Klein Gordon equation that includes negative energies. Particularly, in our formulation both particles and anti- 
particles are modeled on a Hilbert space that includes negative frequencies. We therefore take the point of view that 
our construction is an alternative to the Feynman Stueckelberg interpretation of conventional "in-out" quantum field 
theory. Our approach is similar to standard quantum field theory in the sense that anti-particles are implemented 
as complex conjugates of particles but we do not need to interpret the antiparticles as positive frequency modes that 
"travel back in time". Both the particle field/wave function ip{x) and the antiparticle field/wave function 'fi*{x) are 
propagated with the Wheeler propagator exactly as in the case of real scalar fields. In other words, the time evolution 
of both particles and antiparticles is computed with the same half-advanced half-retarded boundary conditions. 



III. LOCALIZATION 



Locality is generally considered to be a pivotal aspect of quantum field th eory . The Feynman propagator is not a 
local object however. As stressed by for instance Hegerfeldt and Ruijsenaars [20-'25], the time evolution of a positive 
frequency wave function as encoded in the Feynman propagator is such that if the state is well localized initially 
it fills the whole of space an infinitesimal time later. This is one of the issues that has led to quite radical views 
that field theory cannot be considered as the relativistic quantum mechanics of localizable particles, see for instance 
[l3 - [l7j . As explained above however we consider the whole real field to be the wave function. The corresponding time 
evolution is described by the real Klein Gordon equation and the Wheeler propagator and is therefore perfectly causal 
which implies that no infinitely fast spreading of wave functions can occur, see particularly Hence we deduce 
that introducing negative frequency solutions and writing the probability density in terms of the real field itself ([4]) 
cures the non-locality problem associated with the a-causal spreading of wave functions. However, we remind the 
reader that the non-locality of the Feynman propagator is not a problem of standard quantum field theory itself but is 
purely a problem with its "in-out" description, the "in- in" formalism 23- 271 is causal [3. l28l - [30| although not always 



manifestly so, see for instance the Feynman diagram like formulations 31-33|. For applications see e.g. (s^-fs 



A. The Newton- Wigner position operator 



While viewing the real field itself as the wave function eliminates the non-local spreading of wave functions, another 
source of non-locality remains that obstructs a proper first quantized picture. This other non-locality is hidden in 
the localization problem associated with the relativistic quantum mechanics interpretation of quantum field theory. 
The problem is to find a suitable position operator that can be used to compute the average position given a specific 
quantum state. In a famous paper Newton and Wigner [l[ introduced a position operator that is named after the 
authors^. The Newton- Wigner operator is defined on the space of positive energy solutions and therefore inherently 
involves infinitely fast spreading wave packets. It is possible to define analogous position operators on the Hilbert 
space that includes negative energy solutions, see for example 38]. Such operators still possess non-local behavior 
that is inherent to the Newton- Wigner operator though. 

Before discussing our position operator we start by a derivation of the Newton- Wigner position operator for scalar 
particles which differs considerably from the presentation by Newton and Wigner but is very concise and logically 
consistent. We list the well known problems associated with the Newton- Wigner operator. Then we derive a position 
operator that satisfies some of the desirable properties of the Newton- Wigner operator but does not share its problems 
and furthermore provides a covariant extension. 

If one takes the positive frequency Klein Gordon density (jlj seriously as a probability density one is inclined to 
define the position operator through the following expectation value, 

(f) f d^xx (j)X{x)idoMx), (19) 



Pryce was actually the first '45*1 to introduce the Newton- Wigner position operator in the context of photons, but since Newton and 
Wigner gave a systematic derivation from seemingly reasonable axioms their names became associated with the operator. 
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this definition seems obvious but poses problems when going to momentum space. The naive position operator that 
realizes the above expectation value in momentum space is precisely the non-relativistic operator, 

XNR idj:, (20) 
so that the expectation value in momentum space is given by, 



r+ik)id^Mk) ■ (21) 



Because of the non-trivial measure factor (2wfc) ^ in (j2ip the non-relativistic position operator is obviously not 
Hermitian. The Newton- Wigner operator is simply the Hermitian part of the non-relativistic operator, 



^NW ^ i(% + dl) ^ + ^-d~) ^ ^). (22) 



Hence this operator contains a non-trivial function of the momentum, which is at the heart of the problems of the 
Newton- Wigner position operator. First of all this term prevents a simple interpretation of the position operator in 
position space, its representation in position space does not lead to the simple expectation value P^ . Instead the 
operator becomes highly non-local in position space^, 

1 9- 

^NW^{x + --^rf^)- (23) 

As is widely discussed in the literature this non-local behavior makes the interpretation of the Newton- Wigner operator 
as an exact position operator highly questionable. The expression is moreover not covariant and transforms non- 
trivially under Lorentz transformations, so to what extent it can be interpreted as a relativistic operator is not 
clear. We will not discuss particles with spin in depth in this article but we just mention that the Newton- Wigner 
operator does not even exist for massless particles with spin greater then one-half such as photons. Generalizations 
of the Newton- Wigner operator exist that apply to photons but are not physically acceptable since they either have 
non-commuting components or do not transform as three vectors. 

The absence of acceptable relativistic position operators presents us with an awkward situation. Planck and 
Einstein realized that electromagnetic radiation should be emitted in finite energy wave-packets to resolve the ultra 
violet catastrophe of the Rayleigh- Jeans law. The reason why the electromagnetic field should be subdivided into 
quanta is explained by the statistical interpretation of Einstein and Born. The fact that the square of the field 
amplitudes should be interpreted as a probability density enforces the fields to "square integrate" to multiples of 
one, which is a constraint that supplements the "classical" field equations. The unavailability of suitable relativistic 
position operators however implies that there is no probabilistic interpretation for relativistic fields that is similar to 
quantum mechanics. Often this fact is considered to be a mere curiosity but it should be realized that the absence 
of a first quantized Copenhagen like interpretation of relativistic fields removes the naturalness of the partitioning 
of Klein Gordon fields, Dirac fields and Maxwell fields into packets as required by the photoelectric effect and the 
Bose-Einstein, Fermi-Dirac, and Planck distributions respectively. 

In the next section we introduce a position operator that is covariant, local, causal and furthermore commutes 
with the spin generators and is therefore also valid for fields with spin. Such an operator restores the probability 
interpretation of non-relativistic quantum mechanics in a relativistic setting and therefore also restores the naturalness 
of the subdivision of fields into quanta. 

B. The bilinear space-time operator 

We simply define the space-time position operator from the expectation value of its non operator valued components 
in position space. 



d^x x^" (f>{x)do(pH{x). (24) 



Formula II23II is the diflferential analogue of the integral formula presented by Newton and Wigner [l[ . 
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Note that we implement our probability density which includes negative energy modes to maintain causality, the 
computation is quite similar however if one chooses a positive frequency Hilbert space. We use the following off shell 
Fourier transforms, 



d^fc 



(t){k)e 



ikx 



(27r)4 
where we have that 

= -i27rd{k^ + m?)e{k°)f*{k). 



_d4fc' 

W) 



-ik' X 



(25) 



(26) 



Let us start with the derivation of the position operator. The essential step that avoids the problems of the Newton- 
Wigner operator is to represent the components in ()24p by a bilinear derivative in momentum space. 



- / d'x 



d-^fc d'^k' 
(27r)4 (27r)' 



>(fc)0Jf (fc') (fco + k'^) e^('=«"'=o)*(i% - id^)e'^^-^'y\ 



(27) 



Upon evaluation of the on shell delta functions and using partial integration the expectation value can be written as 
follows. 



d^fc 



(27r)32wfe 



/* (fc+ , t)idj:f{k+ ,t) + r (fc- , t)t d j:f{k^ , t) 



which can be written in a covariant manner as, 

d'^k 



i2ny 



rik,t)df:f{k,t) 



/(fc,i) = /(fc)e''=°* 



(28) 



(29) 



Here we stress that in the above expression the integral over the frequency k^ should be performed before taking the 
derivative. Hence, the position operator is given by. 



(30) 



For general wave functions we can see that the average position, or "center of probability" , is given by, 

(f)(i) = (^)(0) + (^3)i. (31) 

Observe that the motion of the wave packet is a solution to Newton's law of motion for particles where an acceleration 
term is absent since the wave packet is free. The initial location of the "center of probability" of the wave packet is 
given by. 



nk+)^^f:f{k+) + r{k^)^^J:f{k^) 



(27r)32a;fe 

and the velocity of the wavepacket is the expectation value of the group velocity. 



d^k dwk 



(|/(A:+)|^-|/(fc_)P). 



(32) 



(33) 



Interestingly, formula ([33]) allows us to understand the group velocity as an operator in the momentum space repre- 
sentation of wave functions. Equation ([5T|) shows that the "center of probability" of general wave packets follow a 
particle like trajectory which should be interpreted in a similar spirit as Ehrenfest's theorem. 

Because we simply define a covariant space-time operator from the expectation value (I24p we also need to promote 
time to a quantum mechanical operator. While conflicting with popular belief, time can be promoted to an oper- 
ator quite easily, provided one takes a Lorentz covariant perspective^. If one chooses to Fourier transform all four 
coordinates to compute the expectation value, time will automatically need to be represented as an operator. 



rf'^ (^^0(fc)0^^(fc') (fco + K) e^^~^-~''y^{^^,, - zafeJe^('=-'=o)' 



(34) 



As operators time and position are on equal footing but the expectation values of these operators are certainly not on the same footing. 
Accordingly, the energy time Heisenberg relation AEAt does not have the same exact status as the momentum position uncertainty 

relation APAx > -^fi. The asymmetry can be traced to the fact that the integration in expectation values is over a spacelike hypersurface 
in a nevertheless covariant manner. 
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If one evaluates the above expression one simply obtains, 

(I) = t, 



(35) 



where t is simply the time coordinate, or equivalently the time parameter, that labels the spatial hypersurfaces. The 
result is simple because we chose our spacelike hypersurfaces to be surfaces of constant coordinate time t. One is 
however free to use different hypersurfaces, then the expectation value of the time operator would be non-trivial and 
give the average coordinate time of the hypersurfaces. Note that this use of the time operator is not directly related 
to time of arrival problems^*^. By partial integration we see that the expectation value of the time operator can be 
written as 



d^fc 
(2^ 



r{k,t)idkj{k,t) ~ r{k)idkj{k) 



(36) 



where we stress that unlike in the case of the position operator one should evaluate the on shell delta functions after 
taking the derivative. We read the above formula as, {i)t = {t}t — {i)t=o such that the time operator is given by 



(37) 



and to = {t)t=o is the value of time at the spatial hypersurface for which t = 0, where the time operator is defined 
such that 



(i)t = t + to = -- 



(2^ 



S{k^+m^) 



r{k,t)tdkof{k,t) 



(38) 



Curiously there is a relative sign if one compares the time operator with the position operator. Combining the above 
results (1501) and ([57]) defines the covariant bilinear space-time position operator. 



-Oh 

2 ^ 



dk 



and the expectation value of the space-time position operator is given by. 



d^fc 
(2^ 



r{k,t)zdkjik,t) 



(39) 



(40) 



Where again we emphasize that, when one computes the position expectation value, one should first perform the 
frequency integral and then apply the derivative while for the time operator the order should be reversed. Furthermore 
it is important that when we write the bilinear derivatives we understand them to act on the wave packets only and 
not on the delta function. Interestingly, we have represented time both by an operator and by a parameter, the 
two properties are not mutually exclusive as is sometimes argued. The expectation value of time as an operator 
simply returns the value of the time parameter that labels the spatial hypersurfaces. The bilinear action avoids the 
problematic non-local and non covariant terms in the Newton- Wigner operator ([221 US]) • 

A similar bilinear position operator was, as far as we know, introduced into the conventional positive frequency 
formalism of quantum field theory by Olkhovsky and Recami, see the for example the recent paper [s^ and the 
references therein. Our operator differs however from that of the aforementioned authors in two essential aspects. 
Firstly, they interpret their bilinear operator as the antisymmetric part of a total non covariant Newton- Wigner like 
operator. The antisymmetric position operator is considered to be a "mean position operator" while the symmetric 
operator describes in a rather opaque manner the width of a relativistic wavepacket. One problem with that point 
of view is that it leads to a peculiar modification of Heisenberg's uncertainty relation. We on the other hand show 
that the width and other characteristics of wave packets can be computed completely according to the standard rules 
of quantum mechanics by the standard deviation and higher order moments. In particular, we show below that our 
interpretation leads to the same form of the Heisenberg uncertainty relations as in non-relativistic quantum mechanics. 

The second aspect where our operator differs is that our time operator is defined as the conjugate of the off shell 
frequency and not as the conjugate of the (on shell) energy. One implication of this is that even for systems with 



Our time operator can also be used to study time of arrival type problems if one uses a different recipe for computing expectation values, 
this is however beyond the scope of the current paper. 
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discrete energy levels time will be represented by a continuous operator. Another consequence is the different sign of 
the time operator and position operators. Moreover, since our time operator is dual to the frequency and we include 
negative frequency modes in our wave functions the exponentiated time operator is a unitary translation operator. 
Time operators that are dual to the energy cannot simply be exponentiated to a unitary operator since the energy 
is bounded from below. This unitarity problem plagues most other proposals for time operators and was discussed 
initially by Pauli p3 |. In other words, not only causality but also unitarity indicates that wave functions that include 
negative frequency modes should be used . 



C. Powers of the space-time operator 



We define the "position-squared operator" from the expectation value, 
1 



(fix x^x^ (j){x)do(pH{x), 



(41) 



where the indices i,j are the spatial components of x^. The position components can similar to before be given a 
bilinear representation in momentum space. 



1 f d^fc 



{2n) 



(42) 



where it should be noted that one should evaluate the ko integral before taking the derivative. Equation (|42|) yields 
an, as far as we know new, real and symmetric bilinear position squared operator, 



(43) 



Using f{k,t) — /(fc)e*'^°' we see that the dynamics of the position squared operator at non zero times behaves as 
follows, 



{x^x^){t) = {x^x^){Q) + {{vlx^) + (vlx^)) t + {vlvD t\ 
where the initial expectation value is given by. 



{x'x^){0) 



d^k 



(27r)32cjfc 

The coefficients of the terms linear in time are of the form, 



+ (+^-), 



1 



d^fc 



^fl(k)^^,J+{k) 



' s ' 4 7 (27r)32a;fc 

and the term quadratic in time is the expectation value of the group velocity operator squared. 



d^k 



(27r)32wfc 



dki dkj 



r+{k)u{k) 



(44) 



(45) 



(46) 



(47) 



With the position operator and the position squared operator at hand we can compute the dynamics of the width of 
a general wave packet. The width in the z direction is for example given by. 



Az^(t) = {P) - {zf+2{{vlz) - {vl){z))t+{{(vlf) - {vlf) 

where the averages are defined with respect to the wave packet ai t — Q which is f{k). 

We define the time squared operator in a similar fashion as the position squared operator. 



(t2 



1 f d'^k 



(27r)3 



5{k' + m')r{k,t) (9fcJ^ + (9fc„)^ f{k,t), 



this can be seen to give, 

(r2) = (t+to)2. 



(48) 



(49) 



(50) 
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where we have that the initial time parameter to is consistent with our time operator, 

^o = -\J + m')rik)^^U{k), (51) 

and we have that, in addition the fact that is simply the square of the above, 

to--^/ ^<5(fc2 + m2)r(fc)(^(g,j2^(g,„y/(fc). (52) 

This follows from our requirement that in position space, 

{t^) ^IJd^'' {t + tof (bix)do(l)Hix) - {t + tof, (53) 

which for t = directly implies ([5^ . 

The time squared and position squared operators can be combined in a covariant operator 



The covariant position operator and its higher moments possess many desirable properties. One important aspect 
is that they make sense as quantum operators since they are Hermitian. The Hermiticity is even manifest as the 
position operator clearly is imaginary and antisymmetric and the position squared operator is real and symmetric. 

D. The Weyl Heisenberg Poincare algebra 

The bilinear position operator is part of a remarkable algebra that embodies both quantum mechanics and special 
relativity without resorting to "second quantization" . The position operator satisfies covariant commutation relations, 

= =0, (55) 



[i^p,] (56) 

This algebra is a relativistic generalization of the Weyl- Heisenberg algebra. To derive the commutation relations test 
functions on the left as well as to the right of the operators are needed since in our setting operators are understood 
to have a bilinear action. The identity generator is related to phase transformations and is a central element of the 
algebra since it commutes with all other generators, 

[^^/] = [p^,/] = [^^.,/]-o. (57) 

The bilinear position operator also gives a differential operator representation of the Poincare algebra. Translations 
in momentum space are encoded in the relations, 

and the standard Poincare algebra is given by, 

[L^^ , Pp] = ?7^p - Vi^pPf, , (59) 



[-^/Jll'l J -^^21^2] Vpit^2 LuiU2 ^flP2 ^^1^12 ^A'11^2 ~^ 'nviV2 -^/Jl/J2' (6^) 

The complete algebra might be named the "Weyl- Heisenberg- Poincare" algebra. It is the algebra of "events" in space- 
time and is the commutative version of Snyder's algebra [40|, which is well known in the field of non-commutative 



geometry, see for instance |41ll42| . In this sense our position operator unifies the coordinate concept of special relativity 
with the position concept of relativistic quantum mechanics. If one uses the Newton- Wigner position operator on 
the other hand one needs to clearly distinguish those two concepts as recognized for example in 'isj. An important 
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lesson that we extract from the biUnear position operator is that the differential operator representation of the Weyl- 
Heisenberg-Poincare algebra should be taken seriously as a defining characteristic of relativistic quantum mechanics 
since it can be consistently realized by Hermitian operators both in momentum as in position space. Moreover, by 
exponentiation of the generators, the Poincare group can be unitarily represented in momentum and position spacc"'^^. 
The standard realization of the generators in position space is 

Lf_iu = Xfj^pi, —ppXfi, — ^ Xfi, pf^ —idf^i. (61) 

The differential representation of the generators in position space (|6ip is well known and their algebra is sometimes 
interpreted in a quantum mechanical sense but one encounters problems in momentum space if one uses the Newton- 
Wigner position operator. It does not satisfy the covariant commutation relations the algebra cannot be realized 
in momentum space. Hence, a complete relativistic quantum interpretation of the Weyl-Heisenberg-Poincare algebra 
cannot be given if one adheres to Newton- Wigner localization. The bilinear position operator does allow us to represent 
the algebra in momentum space however, 

i 

Lf,u = x^p^ - p^x^, x^.^-'dk^, P^^P^ (62) 

Consequently, this realization completes the single particle relativistic quantum mechanical interpretation. We em- 
phasize again that the bilinear position operator should act on wave functions that include negative frequencies, 
otherwise problems with non-locality will remain if one considers amplitudes, and the time operator could not be 
exponentiated to a unitary operator. We also point out that observables do not depend on whether one discusses 
these in position or in momentum space. 



Particles with spin 

At the level of the algebra it is not difficult to introduce spin by adding to the generators of orbital Lorentz 
transformations the spin generators of local Lorentz transformations^^. 

J^iu — L^jj + S^j^. (63) 

If one replaces the orbital Lorentz generators L in (|55] - |60| with the total angular momentum operator J^ji, one obtains 
the correct algebra for spinning relativistic particles. The spin generators seperately generate a Lorentz algebra which 
is formally equivalent to (j60p . and commute with the orbital generators, 

[L^,,Sp„]^{). (64) 

The local nature of the spin generators is stressed by the fact that they commute with translations and coordinates, 

[V,Pp]=0, [V>i'']=0. (65) 

These commutation relations also make it clear that the bilinear position operator is defined independent of the spin 
and mass of a quantum field. This is diametrically o ppo site to for example the Newton- Wigner operator for the Dirac 
field and the Pryce operator for the Maxwell field [43 which depend non-trivially on the spin generators. Pryce's 
operator for the photon is for example closely related to the spinless Newton- Wigner operator 



fc , fc X 



s 



+ (66) 



2-^ 



where in this case cuk — vk^ since photons are massless. The components of the local Lorentz algebra valued vector 
s are the spin one generators of the so(3) sub algebra given by the epsilon symbol {si)jk — Sijk- Correspondingly, the 



As discussed above, even our time operator can be exponentiated to give a unitary operator since it is dual to the frequency, not the 
energy. Since the frequency unlike the energy takes on values over the whole real line, it is not bounded from below, therefore our time 
operator does not suffer from the unitarity problem addressed by Pauli |44|| . 

It is sometimes asserted that spin is an inherently quantum phenomenon. We do not agree with this statement and simply view spin as 
the conserved charge related to active Lorentz transformations in the rest frame, i.e. the local Lorentz transformations. 
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Pryce position operator interacts non-trivially with the spin of the photon. Moreover, the non-trivial spin dependence 
causes the unphysical effect that the Pryce coordinates of a photon do not commute, similar non-commutativity occurs 
for other operators as well, see for instance [H, \^ . 

As emphasized above, the bilinear position operator does not act non-trivially on the local Lorentz spin indices, in 
fact it does not contain any of the other generators of the Lorentz group. It therefore carries no dynamical information 
which implies that the bilinear position operator can be consistently understood as an operator in the Schrddinger 
picture while the Newton- Wigner operator behaves like an operator in the Heisenberg picture. This is exemplified 
by the fact that the commutation relation of the Newton- Wigner position operator with the zero component of the 
momentum has the form of a Heisenberg equation, 

^^v,^^[po,xf'^], (67) 

where v denotes the Newton- Wigner velocity operator. The bilinear three-position operator is time independent on 
the other hand and commutes with the zero component of the momentum as demanded by the Weyl-Heisenberg- 
Poincare algebra. Covariance of the bilinear position operator formalism can therefore be understood as coming 
from a clear separation between kinematics and dynamics. The kinematics is covariantly described by the Weyl- 
Heisenberg-Poincare algebra and the dynamics is described in a covariant manner by the Klein Gordon, Dirac and 
Maxwell equations of motion. 



E. The Heisenberg uncertainty relation in relativistic quantum mechanics 

Since we have constructed a position squared operator as well as a position operator we are able to compute 
the coordinate space width and the momentum width of a wave function as described by the respective standard 
deviations. Using these standard deviations we can derive Heisenberg's uncertainty relations for relativistic quantum 
mechanics in almost the same manner as in non-relativistic quantum mechanics. To make the presentation as clear 
as possible we present Heisenberg's uncertainty relation for covariant relativistic wave functions in 1 -1- 1 dimensions. 
The average position and the average position squared are simply given by 

If ^ If ^ 

{x) — — dx X (f){x)do(l)H{x), {x^) = - dx x^ (l){x)do(f>H{x), (68) 



2 J ' ^ ' " ' ' 2 _ 

Upon inserting the Fourier representation of the field we obtain 

^^^-U (2^^*(^)^^^+(^) + - - -I J (2^^*(^H90V+(P) + (+ - -), (69) 

and the analogous width of the momentum space wave packet is defined by the momentum expectation values, 

= (2^ ^ + '^-^^^''^ ' ^^'^ ^ ' J (2^ ('^+(^^'' + ' (^^^ 

the standard deviations are defined as usual, 

{Axf = {x') - {xf (Ap)2 = - (fc)2. (71) 

The first step to obtain Heisenberg's inequality one uses the Cauchy Schwarz inequality. The inequality holds exactly 
as in standard non-relativistic quantum mechanics, the non-trivial relativistic measure does not lead to alterations, 

(x2)(fc2) > {xkf, (72) 

Also the second step goes through exactly as in non-relativistic quantum mechanics. The expectation value of xp is 
always bigger than the expectation value of its imaginary part —i^[x, k], 

{xkf>l{\\[x,k]\\)\ (73) 

shifting the operators on the left by the average produces Ax'^Ak^ and leaves the commutator invariant. Again 
exactly as for the position operator of non-relativistic quantum mechanics we have the commutation relation [x, k] — i 
and taking the square root on both sides we have that Heisenberg's uncertainty principle holds precisely as in non- 
relativistic quantum mechanics, 

AxAk > i. (74) 
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Equivalence to Schrodinger position operators 

In momentum space one can relate our wave functions /(fc), which are Lorentz scalars, to wave functions that are 
normalized as in non relativistic quantum mechanics, 

4'+{k) ^ ^=J==f{+oJk,k), ^-{k)^^=l==f{-LOu,k). (75) 

With such a redefinition it is not immediately obvious that our position operators are equivalent with the position 
operators of non relativistic quantum mechanics. A short calculation shows however that 

^(2^^/*(±^fc,fc)9fc/(±Wfc,^) = \ {r±{k)du4'±{k)) , (76) 

which follows by antisymmetry of the bilinear derivative, partial integration and the definition of the 9 operator. The 
right hand side of (|76p has the form appropriate for a standard quantum mechanical expectation value. The position 
squared operator is also equivalent to the standard operator of non relativistic quantum mechanics. Up to certain 
boundary terms that vanish at infinity we have that, 

* ^ r{±Uk,k)(fk+B\)f{±u:k,k)=r±{k)dli:±{k). (77) 



2 (27r)2wfe 



So we conclude that the expectation values of our relativistic position operators can be written in the same form as 
in non-relativistic quantum mechanics. Hence it is not surprising that Heisenberg's uncertainty relation is exactly the 
same in the relativistic and non-relativistic cases. Note that we kept the negative frequency contributions explicit 
however. 



IV. LOCALIZED STATES 



An important part of the work by Newton and Wigner [i| is their definition of a so called "localized" state. Their 
"localized" state appears as an eigenstate of the Newton- Wigner position operator. Notice that we used quotation 
marks to indicate that the Newton- Wigner wave functions are in fact only approximately localized. In momentum 
space the Newton- Wigner "localized" wave function is given by. 



f+^iLJk,k) = ^^\/k^+m^. (78) 

Besides being not perfectly localized the Newton- Wigner wave function is also not a proper wave function since it is 
not normalizable. This was already noted by Newton and Wigner and is obvious if one inserts in the Klein Gordon 
normalization formula (|A6|) . The position space form of the wave function makes it clear that it is not perfectly 
localized, 

1 / Tfl \ ^/^ 

f^wit = 0,x)= - if5/4(mr), r = |l.f|l, (79) 

247r2 1 (ij \ r / 

around r — the Newton- Wigner state behaves as -gW. This behavior remains valid in the massless limit, 



1 



j.5/2 • 

5/2 



/+^y(t = 0,f) = -^ - , (80) 
22772 \r J 

Note furthermore that the Newton- Wigner wave function is a positive frequency solution to the Klein Gordon equation. 

In 01 we introduced a wave function for massless scalar particles that unlike the Newton- Wigner wave function 
is real, normalizable and perfectly localized in the sense that in a certain limit the probability density is given by a 
Dirac delta function^'^. For clarity we discuss here the wave packet that is peaked around the origin, both in position 



Our wave function is essentially the real part of the state discussed by G. Kaiser Il9ll . The author interprets his wave function as a 
coherent state, we show in this section however that this state is not a coherent state, it does not saturate the uncertainty relations. A 
similar state also appears in [isll but the authors do not normalize their states in a relativistically invariant manner which allows them 
to use the localization concept of non-relativistic quantum mechanics. However, by using wave functions that are not normalized in a 
covariant manner, they completely ignore the issue of relativistic position operators. Their results are nevertheless qualitatively similar 
to ours since our bilinear operator is the straightforward analogue of the position operator in non-relativistic quantum mechanics. 



14 



and momentum space. We also choose the initial time of the state such that t = to = 0. The wave packet of the state 
in momentum space is given by, 

/+ (fc) = /_(fc)= 4^6 e-^ 11*^11. (81) 
In position space the wave function is defined by, 

6(X) = (4Tre) [ ^-^llkUt-^e)+^k■S ^ g.||fe|l(t+«)+.fc.j\ ,g2) 

J (2^)3||fc|| V / ' 

which gives, with r = \\x\\, 

The Hilbcrt Transform of the wave function is, 

1 / Z6 Z6 \ 

At t = the probability density as defined by dU is, 

P{x) = i0(x)ao0«(x) = ^ l^j^- (85) 

As observed in [2l| the probability density that corresponds to our state is perfectly localized in the limit where e — > 0, 
limP(a;) = (5(3'(f). (86) 

e— >0 

The point we wish to make here is that this sense of localization corresponds precisely to the concept of localization 
as described by our bilinear Hermitian position operators. It is easy to see that since our wave packet in momentum 
space is real and our position operator is antisymmetric, that the average position of the wave function is zero. The 
width of the wave function in position space can be computed by our real and symmetric position squared operator. 
In the z direction we have that^*, 

(27r)'^||fc|| \ / J (2TT)'^\\k\\ 

This integral can be computed for arbitrary positive epsilon by introducing spherical coordinates in the momentum 
space variables, giving (z^) — e^. Since (z) = we have that the width of the wavepacket as defined by its standard 
deviation is simply given by, 

Az = e. (88) 

From this it is clear that the wave packet becomes more and more localized as e approaches zero. So indeed our 
symmetric position squared operator is perfectly compatible with the localization as described by This should 

come as no surprise since our position operators in momentum space simply correspond to the averages of the 
coordinates in the position space expressions. The expectation value (z^) as computed in momentum space (j87p is 
simply equivalent to the position space integral. 



a". 



where P{x) is the probability density given by ([55)) and the integral is solved by going to spherical coordinates. We 
reiterate that the Newton- Wigner localized state does not posses a probability density since it is not normalizable 



Since the state that we discuss is spherically symmetric any other direction would yield the same result. 
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and therefore does not correspond to a perfectly localized wavepacket even if we were to use our bilinear position 
operators. 

The momentum localization properties are described by the standard deviation of the momentum operators. Since 
our probability density in momentum space is even in the momenta the average momentum is zero. The momentum 
squared average can be computed and is given by, 

(kl) = (Anef f ^!^fcfe-2^ll^^ll = ^. (90) 
J (2^)3||fc|| " 2e2 

The width of the wave packet in momentum space as given by the standard deviation is therefore, 

At = -Li. (91) 

The phase space localization properties of the wave function are therefore independent of e, 

AzAfc, - -i=, (92) 

which exceeds the bound set by Heisenberg's uncertainty relation by a factor \/2. The state is therefore reasonably 
localized in phase space but is not optimally localized for any value of e. It is nonetheless perfectly localized in position 
space if e tends to zero and becomes perfectly localized in momentum space in the limit where e tends to infinity. 
For t the probability density is given by. 



4 e^{t'^ +r'^ 

-{t ~ ze)2 + r2]2[_(t + ie)2 + ,^2 



^\^) = ~'\ U ,-,\2 I „212r li- I „-,\2 1 ^212- ^"^1 



The state can of course not be perfectly localized around a; = for arbitrary times by taking a limit. Interestingly 
however the state spreads but remains localized in the sense that is confined to the lightcone, 

lhnP(x) = ^[<5(t-r)+<5(i + r)]. (94) 

Observe that the magnitude of the probability decreases with the area of the spherical surface around the point of 
initial localization. Notice furthermore that the probability density is an expanding spherical wave for t > and a 
contracting spherical wave for t < 0. This indicates that both advanced and retarded solutions to the field equation 
contribute to the probability density since the retarded solution of the field equation is an expanding spherical wave 
valid for t > and the advanced solution is a contracting spherical wave solution for t < 0. This also fits perfectly 
with the half-advanced half-retarded nature of the Wheeler propagator. 

Note that in particular in classical electrodynamics one often only works with the retarded contribution. Without 
further input however such solutions violate conservation of energy rather dramatically, and in our case conservation 
of probability would be violated as well. Our probability density is however perfectly conserved for arbitrary e and 
its limit yields the half-advanced half- retarded result ([M]). 

So far we have only discussed the probability density, it is also illuminating to study the energy density. The total 
energy of our state is most easily calculated in momentum space where it has the simple form of an expectation value 
(lB8l) . for i = we have, 

This value is also valid for arbitrary times by energy conservation. From ([95|l we see that a perfectly localized particle 
seems possible in principle but in practice it would require an infinite amount of energy to do so. Again this agrees 
well with physical intuition based on Heisenberg's uncertainty principle. 

Note that in many respects the wave function that we discuss resembles a Gaussian. It is however difficult to 
compute the relativistic dynamics of a Gaussian wave packet while the time evolution of the wave function discussed 
here can be fully treated by analytical methods. The case of massive fields can also be treated in some detail, we 
chose the massless example however to illustrate our main points related to localization as clearly as possible. 
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V. CONCLUSIONS 



We gave an interpretation of free quantum field tlieory as relativistic quantum meclianics by including negative 
frequency modes in the wave functions. The appropriate single particle Hilbert space for neutral scalar particles can 
be constructed from the probability density, 

P{x)^^cP{x)do<PH{x), (96) 

where ipiiix) denotes the Hilbert transform of 4'{x) with respect to time, where we emphasize that (pnix) is real. 
The probability density ([M)) gives a wave function interpretation to the real Klein Gordon field so the Klein Gordon 
equation becomes the relativistic analogue of the Schrodinger equation. The propagator that replaces the Feynman 
propagator is the Wheeler propagator which is the half-retarded half-advanced propagator, 

(j){x) = ^e{t - t') j d?x'l\w{x - x')dt4{x') + ^e{t' -t) j d'^x'(j}{x")dt'^w{x" - x). (97) 

Since the Wheeler propagator vanishes outside the light cone the dynamics of the real wave function (l){x) is manifestly 
causal. Moreover, since both our propagation law and our probability density contain real variables only implies that 
we are able to formulate quantum mechanics in position space using real variables only. Furthermore, we discussed 
the bilinear position operators, 

^ ^ 9 fc, , ^-^-\(^ g g fe„ + g g ) . (98) 

These operators are covariant, Hermitian, local and well defined for particles with spin such as photons. Therefore they 
are a straightforward extension of the non-relativistic position operators in the Schrodinger picture, unlike operators 
related to the Newton- Wigner position operator. Interestingly, the bilinear nature of our operators makes the reality 
of their expectation values manifest. The antisymmetric operators imply that an imaginary part is taken and the 
symmetric operators imply that a real part is taken. Additionally, the bilinear space-time position operator combines 
quantum mechanics and special relativity in the sense that it allows the covariant Weyl-Heisenberg-Poincare algebra 
to be realized both in momentum and position space. 

Particularly, our space-time position operator explicitly contains a time operator that can be exponentiated to a 
unitary operator. A standard argument due to Pauli that the time operator cannot be unitarily exponentiated is 
rooted in the fact that the energy is bounded from below. This problem is evaded simply because our time operator 
is dual to the frequency and not the energy. Furthermore, since our wave functions also include negative frequencies 
no problems with unitarity occur when our time operator is exponentiated. 

Moreover, we have presented an explicit real solution to the Klein Gordon equation 0g that can be perfectly localized 
in the sense that it has a limit where the width of its probability density becomes zero and the density itself becomes 
a Dirac delta function, 

Iimi0,(a;)ao0f(a;)=<5(3)(f). (99) 

We have also shown that the limit cannot be attained in practice however since the total energy contained in the state 
diverges, which is what one expects from physical intuition based on Heisenberg's uncertainty relations, 

E,= j d^xTooix) = i. (100) 

The width of the probability density in position space is given by e, which can be verified using our Hermitian bilinear 
position operators in momentum space. Hence, the degree of localization of a wave function can be obtained with our 
bilinear position operators and perfectly localized relativistic wave functions can be obtained in the limit of infinite 
energy. This should be contrasted with the "localized state" of Newton and Wigner [l| which is not normalizable and 
therefore does not describe an arbitrarily localized probability density. 

Together with Q we offer a consistent first quantization of free quantum field theory because the covariant dy- 
namics of the wave functions is causal and our position operator and probability density yield a local probabilistic 
interpretation. We expect a similar picture to hold also in perturbative interacting quantum field theory. Pair creation 
and other types of interactions will then play a role but those effects should be treated as small corrections to our 
localization scheme provided the coupling constant is small. Again we contrast this with the Newton- Wigner local- 
ization scheme were particles are not even localizable in free quantum field theory. A consistent local interpretation 
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seems especially useful in modern applications of quantum field theory where localization is important such as cavity 
quantum electrodynamics. 

We are currently investigating the quantitative influence of interactions on the here described concept of localization. 
It would also be interesting to generalize our formalism to curved space-times. 



Appendix A: Conventional field theory as relativistic quantum mechanics 

The relativistic analogue of the non-relativistic Schrodinger wavefunction is the real Klein Gordon field since both 
describe particles without spin and electric charge. The real Klein Gordon field contains both negative and positive 
energy modes however and satisfies the second order Klein Gordon equation which make an interpretation of the real 
field as a wave function difficult. The positive and negative frequency solutions separately satisfy first order equations 
however, 



ihdt4>±ix) = ±TOc^y 1 - 92 0^(2.)^ (Al) 

which have the form of a Schrodinger equation. This is especially clear from the non-relativistic approximation where 
the free particle Schrodinger equation is recovered exactly, with the appropriate relativistic corrections, 

±ihdo (l,±ix) = mc^(b± - ^d'^± - l-^dU± - - (A2) 

2m 8 m'^c^ 

Since the positive frequency components of real fields are complex valued one can interpret these as wave functions by 
using a relativistic analogue of the "wave function squared rule" of Max Born^^ in non-relativistic quantum mechanics 
[4a £0]. 

P{x) = (j)*^{x)ido(j>+{x). 

whereas the Schrodinger form of the Born rule is simply, Psch{x) = ip* {x)'ip{x). The Schrodinger probability density, 
and therefore also the wave function, is not Lorentz covariant which is clear from the normalization condition, 

d^xi^*{x)il;{x) = 1. (A3) 

The integration element is not a Lorentz scalar but a four vector valued three form^®, therefore the Schrodinger 
probability density is not a Lorentz scalar, in fact it should be a vector valued current. The covariant nature of the 
probability density can be made manifest if normalize wave functions by a boundary integral, 

d^xd''i(t>lix)idMx))= J d''x^r+{x)td^^+{x)^l, (A4) 

which with a standard choice of hypersurface and coordinates this leads to the probability density ([T|). Note that 
(/)_i_(a;), contrary to ipi^)-! is a Lorentz scalar. The relativistically invariant momentum space wave functions are 
obtained by the standard on shell Fourier transform, 

M^) = I + m'nk")f{k)e^'^, (A5) 

such that similar to the wave functions (j)+{x) the wave packets /(fc) are Lorentz scalars. The normalization condition 
on the wave packets f{k) is simply the Fourier transform of (|A4[) in standard coordinates, 

^6{k' + m^)0{ko)nk)f{k) = J k) = 1, (A6) 



As acknowledged by Born himself in his Nobel lecture [50(. Einstein was actually the first to attribute a statistical interpretation to the 
amplitude squared in the context of photons. It is amusing however that to this day the probabilistic interpretation of photons is far 
from clear, see for instance [sill . 

^® The integration element is in fact (Px^ = ^dx'^ Adx^ Adx" e which in standard choice of coordinates, adapted to a spatial hyper- 
surface, reduces to d^x = dx^dx^dx^, see for instance the discussion on p. 134 of Ryder |52|1 . 
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where ujk — Vfc^ + m? and /+(wfc,fc) is a Lorentz invariant positive frequency wavepacket. The positive frequency 
contribution k) is a wave function in the sense that it gives a probabihty density by the following Born rule, 

P{uJk,k,t) = ^^-l^Jl{u;k,k,t)f+{LUk,k,t). (A7) 

For a text book discussion on wave packets in quantum field theory consult for example section 3.1. of [53l |. section 
3.4. of [131 or section 4.5. of [s^. To summarize one might say that in the quantum mechanical interpretation of 
conventional quantum field theory quantization is essentially the extraction of the positive frequency components from 
the fields and enforcing the normalization condition (|A6p as required by a probabilistic interpretation. Note that this 
can be done without promoting the fields to operators. So from this point of view the operator language merely 
is a convenient tool to organize multi-particle wave functions but not essential for quantization. We used the term 
quantization in the restricted sense of providing the space of solutions to the equations of motion with an appropriate 
Hilbert space structure, see (jA4p and (jA6l) . associated to the positive frequency part of the fields. 

We stress however the well known fact that the above discussed quantum mechanical picture of standard quantum 
field theory is not physical however as issues with causality and non- locality are present. Because of these problems 
the first quantized interpretation is most often abandoned. One should realize that nonetheless it is the interpretation 
used in textbook methods to convert scattering amplitudes to scattering probabilities and cross sections. The modern 
interpretation of quantum field theory is not the one that is used in scattering problems. Most often quantum field 
theory is given a functional interpretation, the field is not considered to be a wave function but an infinite set of 
oscillator degrees of freedom. It is the purpose of this paper to show that the field can be interpreted as a relativistic 
wave function provided one includes its negative frequency modes and uses a suitable position operator. 



a. Time evolution violates Einstein causality 



One of the problems with interpreting quantum field theory as a first quantized theory is that the time evolution of 
positive frequency solutions to the Klein Gordon equation, which are meant to be interpreted as wave functions, do 
not respect Einstein causality 0, [20l [23 | . This can be seen by realizing that the Feynman propagator from standard 
perturbation theory governs the time ordered evolution of the positive energy modes, 

e{t - t')<f)+{x) = i J d^x'Apix- x')dt'(l>+{x'), (A8a) 

eit' ~ t)(j)l{x) = i [ d^x'(t>Xix')dt'AF{x-x'). (A8b) 



Formula (jA8a| shows that the Feynman propagator is the retarded propagator for positive frequency fields and (|A8bp 
shows that it is the advanced propagator for the complex conjugate positive frequency field. Notoriously, the Feynman 
propagators possesses contributions outside of the light-cone which directly implies non-causal spreading of a wave 
function 0+ by the propagation law (|A8|) , for an explicit example see [3] . The contributions that violate causality are 
clear from the explicit form of the Feynman propagator in position space, 

Af{x) = e{-x^) ™ (m^/^{x^ +ie)) + ie{x^) ]^ TC^ (m^/x^ + ie) , (A9) 

8^^^y-{x^ + ie) ^ ' At:'^ V x'^ + ie ^ ' 

where the second term is nonzero outside of the lightcone, which follows from the theta function dix?') ~ 9{—c^t'^ + x'^) 
which enforces spacelike contributions only. 



Appendix B: Causality from negative energies 

Perhaps less often mentioned then is deserved but emphasized in our recent work i2] is the fact that the real part 
of the Feynman propagator is perfectly causal and vanishes outside of the Lightcone. We call this propagator the 
Wheeler propagator following ^56] due to its implicit appearance in the Wheeler Feynman absorber interpretation of 
classical electromagnetism [53| . A more modern application of the Wheeler propagator is the perturbative description 
of the classical gravitational problem of pointlike particles [ssl . [59| . It is also known as the principal part propagator 
since it's momentum space definition involves a principal part The operational significance of the Wheeler 
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propagator as stressed in [5] is that it governs the time ordered time evolution of the whole real field and not just its 
positive frequency part, the Wheeler propagator applied to the field gives 

ie(i-<Xx) = [ d^x'Awix-x')dt'<Pix'). (Bl) 



While the Feynman propagator is related to the a-causal positive frequency Wightman function, the Wheeler propa- 
gator is related to the causal commutator function^^ by time ordering, and is therefore also related to the "classical" 
advanced and retarded propagators, 

Aw{x-x') = h{t-t')Acix~x') + ^0it' ~t)Acix' ~x), = ^A^ + ^A^, (B2) 

where the commutator can be given the integral representation, 

Acix + ™^)^(fco)e^'=(^-^'^ , (B3) 

and e(fco) is the sign function which written in terms of Heaviside functions is 6{ka)—9{—k{)).T}ie integral representation 
of the Wheeler propagator is. 



P 



(B4) 



(27r)4 B + m2 

where P denotes the principal value. The explicit form of the Wheeler propagator which follows from (IB2[) and (jB3 



Avv(x) = ^6{-x^) + e{-x^) ^ "l_ Ji (mv^) . (B5) 
27r Anyf-x^ V / 

From this expression it is manifestly clear that the Wheeler function cannot propagate fields outside the light-cone. 
Note that (jBSP is the real part of the regulated expression of the Feynman propagator (IA9p . and (jB4p is the real part 
of the standard integral representation for the Feynman propagator, in both cases the regulator e is taken to zero. 

Also note that the Wheeler propagator, not the Feynman propagator, is the natural Green's function of a real scalar 
field. This is rooted in the fact that Wheeler propagator is a real and symmetric solution to the inhomogeneous Klein 
Gordon equation. 



1. The in-out boundary conditions implied by the Wheeler propagator 

A perhaps more intuitive understanding of the half-advanced, half-retarded nature of the Wheeler propagator can 
be gained from the following form of the propagation law, 

(t,{x) = ]^9{t - t') j d?x'Aw{x - x')dt4{x') + ]^e{t' ^t) j d^x'(j){x")dt'Aw{x" - x), (B6) 

which is equivalent to (|Bip . The physical interpretation of this equation is as follows, half of the field (t>{x) at time t is 
specified by initial conditions and the other half of the field is specified by "final conditions" . The initial condition is 
the specification of a real field (t){x') at a time t' earlier than t and the final condition is the specification of a field 4>{x") 
at a time t' later than t. The use of initial and final conditions on the field implies the use of two Dirichlet boundary 
conditions as opposed to the half-Dirichlet, half-Neumann Cauchy boundary conditions which are standard when 
viewing the Klein Gordon field as a "classical" field. In the context of quantum cosmology the Dirichlet boundary 
conditions are referred to as "in-out" boundary conditions and the Cauchy boundary conditions are referred to as 
"in-in" boundary conditions, see for instance |27| . 

The Feynman propagator also corresponds to a choice of "in-out" boundary conditions. Unlike our "in-out" bound- 
ary conditions which are compatible with the Wheeler propagator, the Feynman Stucckelberg boundary conditions 



^"^ The commutator two-point function is also referred to as the spectral propagator *23*| in the Schwinger Keldysch, closed time path or 
in-in formalism popular in cosmological applications [24l427l| . 
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are complex even for real fields. The initial conditions in the standard Feynman Stueckelberg formalism are not 
specified by the real field — (t>+{x) + 4>-(x) but by the positive frequency solutions 4>+{x) only and the final 
conditions are given by (x) . Equivalently, by the identity for real fields 4>- (x) = ip*^ {x) one can add the negative 
frequency solutions (j)-{x) to the final conditions and their complex conjugates 01 (x) to the initial conditions. This 
corresponds to Feynman's frequently misinterpreted statement that negative frequency solutions "travel back time" . 
Note that with our boundary conditions there is no difference in the "direction of propagation" between positive and 
negative frequency fields. They only appear in their real combination and therefore both are given half-advanced, 
half-retarded boundary conditions. 



2. Negative frequency modes have positive energy 

Historically, negative frequency modes were thought to lead stability problems due to the misconception that 
negative frequency modes have negative energy. The lore was that negative frequency wave functions would not allow 
a true ground state of lowest energy so that the vacuum would decay indefinitely to lower energy states. Dirac's sea 
of occupied negative energy states is a famous attempt to resolve this problem for Fermions. In the modern Feynman 
and Stueckelberg approach both Bosonic as Fermionic negative frequency states are rewritten as positive frequency 
states with advanced boundary conditions. 

We do not agree with the negative sentiments towards negative frequency modes. As stressed in Q the essential 
point is that the energy contained in a field is not given by the energy that appears in the dispersion relation, the 
true energy contained in a wave function is given by the energy momentum tensor, 

j d^x°no = ^J Sx{do^dad? + dd^-d^ + fi^cb^), (B7) 

-lj-0^6{e + m')\ko\\m\\ (B8) 

Note the absolute value operation on fco, it directly implies that the negative frequency modes do not have negative 
energy. We therefore prefer to use the terminology "negative frequency" solutions as opposed to "negative energy" 
solutions^®. Furthermore, notice that the on-shell formula for the energy (jBSp has the form of a quantum mechanical 
expectation value. 



3. Negative frequency modes and the probability density 

In this section we summarize our discussion in [2*1 which was inspired by Halliwell and Ortiz and define a 
probability density on the space of real fields. As written in (|B7p the energy of a scalar field can be interpreted as a 
quantum expectation value. Hence, the energy momentum tensor suggests the following normalization condition for 
the single particle wave function^^, 

I J 0'5(fc2+-')l/(fc)l' = l> (B9) 

which is almost the same as the standard normalization of the Klein Gordon field (|A6|) except that no Heaviside 
function is employed to restrict the integral to positive frequencies. Consequently, a Lorentz invariant version of 
Born's rule for the probability density follows that includes negative frequency wave modes, 

P(fc) = i 1 dkoSik' + m2)|/(fc)p = Ij^^^^ (l/+(fc)P + \f-ikr) . (BIO) 

In position space this gives a negative frequency completion of the "Klein Gordon Born-rule" ((T|) that despite the 
presence of negative frequencies is positive definite, 

Pi^) = 2 ['l^*+i^)'^9o(j)+{x) - (f>*_{x)ido(j>~{x) j . 



In the literature there do appear attempts to add modes whose energy is considered to be negative [6214641 . We do not follow such an 
approach however. 

This form of the innerproduct is discussed by Halliwell and Ortiz in Q who were inspired by work of Henneaux and Teitelboim Q . The 
innerproduct also appears in a series of papers by P. Morgan where a classical stochastic interpretation is given, see e.g. |65ll . 
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This formula is merely a rewritten form of the Klein Gordon form of the probability density since 0_ = 0^ for real 
fields. Hence, adding negative frequencies does not change the probabilities of standard quantum field theory. Inspired 
by we can write the probability density equivalently as, 

P{x) = i J d^x'(f>ix)dtAH{x,x')dt'Hx') 

where the Hadamard function, sometimes referred to as statistical "propagator" and related to the anti-commutator, 
is given by, 

Ah{x -^')^-J ■0j3^i^' + , (Bll) 

where we introduced an overall minus sign to avoid a minus sign in the definition of the probability density ([3]). Note 
that the Hadamard form of the probability density ([3]) is equivalent to the Klein Gordon density ([1]). The essential 
difference between the two representations is however that in the Hadamard probability the real fields themselves 
are relevant and not the positive or negative frequency contributions separately. This shows that the causal Wheeler 
propagator can be used instead of the a-causal Feynman propagator to compute the dynamics of wave functions and 
probabilities. Notice that the Hadamard propagator and its time derivative only appear with equal time arguments, 
therefore it does not influence the causal propagation of the probability density. The Hadamard form of the probability 
density can directly be compared to the Klein Gordon probability density in the following form, 

f ^ ^ 

P{x)^ d^x'(l>{x)doA+{x-x')dQ'<p{x') , (B12) 



This expression reduces to the form familiar from the position space probability density that we introduced at the 
beginning of this paper ^ . To establish the equivalence of (|B12I) and ^ one uses the fact that the positive frequency 
Wightman function is a propagator for positive frequency fields, 

= J d^x'A+{x-x')do'<l)+ix'), (B13) 

We emphasize that in standard quantum field theory the Wightman, or including time ordering the Feynman, propa- 
gator plays a dual role, it both gives the time evolution of the wave function by means of its propagation law (jB13|) and 
it is used to define the innerproduct (|B12[) . If one chooses however to work with real fields as wave functions, the time 
evolution and the probability density are determined by two different two-point functions. The Hadamard function 
defines the probability density ^ and the time evolution is governed by the commutator two point function, or its 
time ordered counterpart the Wheeler propagator (jBl[) . does not satisfy a law similar to (|B13|) . The commutator two 
point function on the other hand is the appropriate (non time ordered) propagator for real fields. 
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